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Abstract. In this work, we have applied the self-compressed stabilized jellium model 
to predict the equilibrium properties of isolated thin Al, Na, and Cs slabs. To make 
a direct correspondence to atomic slabs, we have considered only those L values that 
correspond to n-layered atomic slabs with 2 < n < 20, for surface indices (100), 
(110), and (111). The calculations are based on the density functional theory and self- 
consistent solution of the Kohn-Sham equations in the local density approximation. 
Our results show that firstly, the quantum size effects are significant for slabs with 
sizes smaller or near to the Fermi wavelength of the valence electrons Ap , and secondly, 
some slabs expand while others contract with respect to the bulk spacings. Based on 
the results, we propose a criterion for realization of significant quantum size effects 
that lead to expansion of some thin slabs. For more justification of the criterion, we 
have tested on Li slabs for 2 < n < 6. We have compared our Al results with those 
obtained from using all-electron or pseudo-potential first principles calculations. This 
comparison shows excellent agreements for Al(lOO) work functions, and qualitatively 
good agreements for the other work functions and surface energies. These agreements 
justify the way we have used the self-compressed stabilized jellium model for the correct 
description of the properties of the simple-metal slab systems. On the other hand, 
our results for the work functions and surface energies of large-n slabs are in good 
agreement with those obtained from applying the stabilized jellium model for semi- 
infinite systems. Moreover, we have performed the slab calculations in the presence 
of surface corrugation for a selected Al slabs and have shown that the results are 
worsened. 
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1. Introduction 

A thin slab is a system composed of a few atomic layers in such a way that it is finite in 
one direction (here we take as z axis), and infinite in the other two directions of x and 
y. The finiteness of the size in the z direction gives rise to the so-called quantum-size 
effects (QSE)P ElEllllElEllTlElinillOllIIlITS]. The QSEs are significant when the 
size, L, of the system in one direction is comparable to the Fermi wavelength, Ap, of the 
electron in that direction and become less important in the limit L > 2Af. 

To study the electronic structure and mechanical equilibrium properties of the 
simple metal slabs of Al, Na, and Cs, we have used the stabilized jellium model 
(SJM)[T3]. It has been shown that the SJM, in which the discrete ions of atoms are 
replaced by a uniform positive charge density, is a simple and realistic model to describe 
the bulk and surface properties of simple metals. Moreover, it has been shown that 
the self-compressed version of the SJM (SC-SJM) is appropriate for finite systems in 
which the surface effects causes the jellium background density differ from that of the 
bulk [T^ [T5l [T6] . In the prior work by Sarria et. al.[9l, the SC-SJM has been applied to 
slabs with different sizes of L. For each fixed L value, they have determined the r* value 
for the background that minimized the total energy per particle [E/N), and was always 
smaller than that of the bulk. That procedure was repeated for different values of L 
parameter. They have shown that for certain values of L which were integer multiples 
of Af/2, the self- compression effects are more pronounced. However, since in their work 
the self-compression procedure was performed for constant L values, any change of the 
rg value was equivalent to change in the number of atoms and thereby, the slab systems 
they studied, would not correspond to realistic isolated slabs. In order to apply the 
SC-SJM for isolated slabs, one has to consider as a fixed parameter and let the value 
of L relax to its equilibrium value For fixed A^, the relaxation of L relative to the 
corresponding bulk value leads to a change in the value. We denote by rj the value 
corresponding to the L\ and will show in the following that under certain conditions it 
can assume values larger than that of the bulk. 

In this work, using the SC-SJM, we have calculated the equilibrium properties of 
isolated thin slabs of simple metals of Al, Na, and Cs. To make a direct correspondence 
to atomic slabs, we have considered only those L values that correspond to n-layered 
atomic slabs with 2 < n < 20 for surface indices (100), (110), and (111). The 
calculations are based on the density functional theory (DFT)[T7j and solution of the 
Kohn-Sham (KS) equations |18j in the local density approximation (LDA)[T8] for the 
exchange- correlation (XC) functional. In our calculations, the slabs are taken as isolated 
systems which can undergo relaxations in the z direction. Our results show that, in 
contrast to the r* values (obtained for fictitious constant- volume slab system), for some 
slabs the r| (obtained for realistic constant- A^ slab system) assumes values that are larger 
than that of the bulk. That is, in some cases the relaxations are realized as expansions. 
The expansion behavior which is predicted in this work, have been also predicted by 
the other first-principles calculations [6[ HOl HH |T2]. The results show that the amount 
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of relaxations decrease for large-n slabs. Based on the results, we explore a criterion for 
realization of significant quantum size effects that lead to expansion of some thin slabs. 
For more justification of the criterion, we have tested on Li slabs for 2 < n < 6 and (in 
contrast to the results of Ref. [9i|) have obtained some expansions. In section |2] we have 
explained the calculational details and in section [3l we have discussed on the results of 
our calculations and those obtained using first-principles methods [21 [TOl [8], [7] as well as 
those obtained from SJM calculations for semi-infinite jellium systems[T9l |20]. At the 
end, in Appendix, we have presented the detailed derivation of the relation between r* 
and rj. 



2. Calculational details 



In the SC-SJM, the energy of a system with electronic charge distribution n(r) and 
background charge n+(r) is given by [HI El]: 

^sjm['^,^+] = EjM[n,n+] + [euirs) +WRirs,rc)] j drn+{v) (1) 
+ {5v)Y^s{rs,rc) I dr e(r)[n(r) - r2+(r)]. 



where, i?jM is the energy in the jellium model (JM)[22j, = — 92;/10ro and wr = 27rnr^ 
are the Madelung energy and the average of the repulsive part of the pseudopotential [23] 
over the Wigner-Seitz (WS) cell of radius tq = z^^^rg, respectively, z and (5f)ws = 
(3r^/2r^ — 3z/10ro) are the valency of the atom and the average of the difference 
potential over the WS cell, respectively [13]. n+(r) = n9(r) is the jellium density in 
which n = 3/47rrf and G(r) takes the value of unity inside the jellium background and 
zero outside. Tc is the core radius of the pseudopotential and is determined in such a 
way that the bulk system becomes mechanically stabilized at the experimental value. 
(All equations throughout this paper are expressed in hartree atomic units.) For a finite 
system (with fixed number of particles) at mechanical equilibrium, the value of the 
background is different from that of the bulk and assumes the value for which 

d(E/N) 

^ / ^ -0, (2) 



drs 

where is the number of electrons in the system. 

In the SC-SJM, a slab of thickness L has a full translational symmetry in the x and 
y directions, and therefore, the physical quantities depend on the spacial z coordinate. 
Moreover, the KS equation reduces to a one-dimensional equation, 

V ~2 + ) i'niz) = Sn'lpniz), (3) 

where 

Veff{z) = <f){z) + v,,{z) + {5v)^sQ{L/2 - \z\). (4) 
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Table 1. Bulk values, lattice constants, and interlayer spacings in atomic units. 
The values in parentheses are in units of corresponding bulk Fermi wavelengths, Ap- 



element structure 



a 



dm 



Al 
Li 
Na 
Cs 



fee 
bcc 
bcc 
bcc 



2.07 7.64 (1.13) 3.82 (0.56) 2.70 (0.40) 4.41 (0.65) 

3.28 6.60 (0.62) 3.30 (0.31) 4.66 (0.44) 1.90 (0.18) 

3.99 8.10 (0.62) 4.05 (0.31) 5.73 (0.44) 2.34 (0.18) 

5.63 11.43 (0.62) 5.72 (0.31) 8.08 (0.44) 3.30 (0.18) 



(j){z) is the electrostatic potential energy, 

(j){z) =A7T [ dz' {z - z')[n+{z') - n{z')] 



(5) 



— oo 



satisfying 0(±oo) = 0, and the density is given by^: 



(6) 



with the Fermi energy determined by the charge neutrality condition. 

To solve the KS equation for a slab of width L, we expand the KS orbitals in terms 
of the eigenfunctions of a rectangular infinite potential well[3] of width D L + 4Af, 
which allows a 2Af vacuum at each side of the slab. Our testings shows that this choice 
of the width D ensures the independency of the results to the presence of the walls. 
The KS orbitals are expanded in terms of 150 lowest eigenfunctions of the rectangular 
box of width D. This number of basis ensures the reliability of energies to, at least, six 
significant figures. 

For an X{hkl) slab with X=A1, Na, Cs and (/iA;Z) = (100), (110), (111), composed 
of n atomic layers, the size is given by L = ndhki, with dhki being the interlayer spacing. 
For a unit cell of riat atoms with valency z, the lattice constant a is determined from the 
electronic value by a = (Anzriat/'iY^'^rs. In Tabled] we have presented the numerical 
bulk values of the electronic density parameter r^, lattice constant a and the interlayer 
spacings d^ki for each element. It should be mentioned that since Af = 27r(4/97r)i/3r., also 
linearly scales with r^, the lattice constant in units of Fermi wavelength depends only on 
the valency of the atoms (z) and the crystal structure (riat), i-e., a = {SzHat/SirY^'^XF. 
To proceed with SC-SJM calculations for an isolated slab, we must use a reasonable 
relation between the and L values. Since there is no relaxations in the x and y 
directions, the surface area of the slab remains constant and the change in is solely 
due to the change in L. Here, we assume the simplest relation satisfying the constant 
surface area and constant number of atoms. 



where the pairs {L\ rj) and (L, r^) correspond to the equilibrium and bulk states of 
the slab, respectively. Using this simple relation in our SC-SJM calculations, we predict 




(7) 
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Figure 1. Equilibrium values, in atomic units, for n-layered Al slabs with (111), 
(100), and (110) surfaces. The dotted lines correspond to the bulk value (2.07). 



that in the Na and Cs cases some slabs expand while, all Al slabs contract. Finally, to 
include the surface corrugation, we have used a face- dependent relation^Sjfor {6v), 

2" 



(8) 



d 
To. 

in which d is the interlayer spacing for a given surface. For d and tq, we assume their 
corresponding bulk values. Our calculations show that the self-consistent variations of 
d and tq do not lead to any equilibrium state. Moreover, eventhough the second term in 
the right hand side of Eq. ([8]) makes the work functions face-dependent, our results for 
Al(llO) slabs with n=9, 11, 13, and 15 show that the results are worsened with respect 
to those obtained using a fiat surface. 

3. Results and discussion 

Solving the self-consistent KS equations in the SC-SJM, we have calculated the 
equilibrium sizes, L\ of n-layered slabs (2 < n < 20) for Al, Na, and Cs. For each 
element, the ?7,-layered slab calculations has been repeated for the different (100), (110), 
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Figure 2. Same as in Fig. [T] for Na, witli bulk value 3.99. 



and (111) surface indices. Here, it is assumed that the interlayer spacing is the only 
parameter that differs for differet surface indices. However, as we discuss later in this 
section, by introducing the surface corrugation, {6v) becomes also face-dependent which 
leads to different surface properties for slabs of different indices. The advantages of using 
dhki is that we consider only those L values that correspond to real atomic slabs and 
therefore, for each structure there is a one-to-one correspondence to atomic slabs. For 
each n-layered slab, the KS equations are solved self-consistently for different L values 
and thereby, the value which minimizes the total energy per particle, is determined. 

In Fig. [H we have plotted the r] for Al(lll), Al(lOO), and Al(llO) as function 
of the number of atomic layers, n. (Throughout the paper, all sub-figures from top to 
bottom, are ordered according to a decreasing bulk interlayer spacing.) The results show 
that rj < 2.07 for all Al slabs. That is, the model predicts all Al slabs are contracted. 
Comparison of the three sub-figures in Fig. [1] and taking into account the different 
interlayer spacings from Table [H we see that the smaller the size of the slabs the larger 
is the contraction rate. Moreover, since the sizes of the 2-layered Al slabs are larger 
than Ap, we do not observe any significant oscillations. 

In Fig. [21 we have plotted the r] results for Na, which shows significant deviations 
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Figure 3. Same as in Fig. [T] for Cs, with bulk value 5.63. 



from the bulk value up to n=13 for Na(lll), up to n=8 for Na(lOO), and up to 
n=5 for Na(llO). These sizes correspond to slab thicknesses of ~ 2Af. For sizes above 
these thresholds the oscillations become negligible. That is, the QSEs are significant 
for L < 2Af- On the other hand, we observe that the oscillations become sharper and 
significant as we go from top subfigure to the bottom. Since from top to bottom the 
interlayer spacing decreases, we are led to the conclusion that if the interlayer spacing, d, 
happens to be so small or the Fermi wavelength happens to be so large that the condition 
2X-p/d ^ 1 is met, then there would be a chance that the QSE be realized as consecutive 
expansions and contractions. However, since this condition is not met in the Al case 
{d ~ 0.5Af), we do not observe any significant oscillations. This criterion along with the 
first-principles or experimental results enables one to estimate an effective free-electron 
value for a given slab. For example, as is shown in the following (Fig. Hj), according 
to first-principles calculations, the Al(lOO) and Al(lll) slabs undergo expansions |10j. 
This means that, using the criterion, the actual Fermi wavelengths along those directions 
should be larger than the bulk free-electron value or the actual electron density in the 
mid of two adjacent atomic layers should be much smaller than that of the bulk. 

The r\ results for Cs slabs are shown in Fig. [3l As is observed, the shapes are 
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Figure 4. Relative total relaxations of Al slabs at their equilibrium states for (111), 
(100), and (110) surfaces. The open circles, sohd circles, and upside-down triangles 
correspond to the results of present work, FPLAPW 10 , and PP-LDA|2], respectively. 



similar to the corresponding shapes of Na slabs. This fact is explained by noticing that 
in Cs the interlayer spacings, in units of Ap, are the same as those in Na (see Tabled]). 
On the other hand, since the density parameter, r^, of the electrons in Cs is larger than 
that in Na, the oscillations in Cs case have smaller amplitudes relative to those in Na. 
According to Eq. ([71), the variations in L and are related by 



V L ; - 

which can be used to determine the total thickness relaxations of the slabs. In Fig. |H we 
have compared the results of total relaxations of the Al slabs with those of first-principles 
calculations [TOl [2]. In the top sub-figure, we have compared the results on Al(lll) 
with those obtained [10] from full-potential linearized augmented plane wave method 
(FPLAPW) in the generalized gradient approximation (CCA) [23]. The FPLAPW 
results show expansions up to n=17 whereas our results show contractions. As we 
have argued, this could be as a result of larger actual electron Fermi wavelength in 
that direction compared to that of the bulk. Similarly, in the middle sub-figure, the 
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Figure 5. Same as in Fig. [4] for Li. 



FPLAPW results show expansions, while our results show contractions. However, the 
(110) slabs in the bottom subfigure show contractions as predicted by FPLAPW and 
LDA pseudo-potential (PP) calculations [2]. In this case, the average behavior of the 
first-principles results agrees well with our results, and we can conclude that Ap along 
this direction is not much different from that in the free-electron model. 

The results of total relaxations for Na and Cs behave similar to their corresponding 
rj plots and are not presented. To close the discussion on the thickness relaxations, 
we have presented the relaxations of Li slabs for n=2, 3, 4, 5, 6 in Fig. O The results 
show expansions at n=4, 6 for Li(llO); n=2, 3, 4 for Li(lOO); and n=2, 5 for Li(lll), 
respectively. These expansions are verifications of our criterion for observing significant 
QSEs. In addition, by Eq. (Q, at these n values, the rjs are greater than that of the 
bulk while the r* values, obtained for Li in Ref. were always smaller than that of 
the bulk. 

The work function of a given metal surface is the minimum energy needed to remove 
an electron from inside of the metal to a point outside across that surface. In the 
slab system, it is equal to the absolute value of the Fermi energy. In Fig. [6l we have 
plotted the work function results of Al slabs and compared with those of FPLAPWfTO]. 
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Figure 6. Work functions of Al slabs at their equilibrium states for (111), (100), 
and (110) surfaces as functions of number of layers. The open circles, solid circles, 
open triangles, asteriscs, and upside-down open triangles correspond to the results of 
present work, FPLAPW[inj, PP-LDA[S], PP-LDA[7], and PP-LDAg], respectively. 
The dotted Hues correspond to the SJM results for the semi- infinite system [l9t [20]. 



However, for completeness we have also included the pseudo-potential calculation results 
in the LDA (PP-LDA) [21 [HI E] as well as the SJM results for semi- infinite system [HI [20]. 
In order to explain our results, as is observed from FPLAPW results, in Al(lOO), there 
is an excellent agreement in the values and oscillation behavior for n > 4. On the other 
hand, the agreement is poor for n=2, 3. In other words, the bending of the FPLAPW 
curve is a little bit slower than that in our results. Inspecting the top and bottom 
sub-figures also shows that our results and those of FPLAPW have similar oscillatory 
structures. The relative vertical shifts at asymptotic large-n regions can be related to 
the differences in the values of the surface dipole moments (SDM). The SDMs in our 
model calculations are higher than those in the FPLAPW or PP calculations. It is 
because, in the SC-SJM the background density has the same constant value over the 
slab width, whereas in the FPLAPW or PP the discrete ions of each plane are allowed to 
find their own equilibrium positions. It means that, if we had taken an inhomogeneous 
jellium background to simulate the discrete ionic planes, then the equilibrium densities 
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of background sub-slabs would be different. This behavior will somehow decrease the 
separation between the positive and negative charges and thereby decreases the SDM. 
In the extreme case of complete deformable background, the ultimate jellium mo del [25]. 
there is no charge separations and therefore, the SDM vanishes. On the other hand, we 
should not forget that in the FPLAPW calculations the GGA has been used which gives 
results different from that in the LDA for highly inhomogeneous electron densities at 
surfaces. This fact is shown by presenting the PP-LDA results in the middle and bottom 
sub-figures. The PP results lie above those of the FPLAPW with the same structure. 
Returning back to the oscillatory structure, the local Fermi wavelengths of the electrons 
in the FPLAPW are different from those of SC-SJM. Under the condition {Xp/d) > 1, 
a small decrease in Ap (other parameters kept fixed) results in a decrease of the value 
of the ratio, and in turn, an increase in the bending rate. It is the case in our Al(lOO) 
results. Here, slabs with n=2, 3 have been contracted while those in the FPLAPW are 
expanded (See middle sub-figure in Fig. H]). It should be mentioned that for small n's 
no bulk region exists for the slab (unless Xp/d <^ 1), and expansion or contractions will 
have different impacts on the surface properties. However, for large n's, the effects of 
relaxations in the bulk can be canceled out by the relaxations of the neighboring layers 
and therefore, the SDM does not change. In other words, the relaxation of the layers 
which lie inside the major peak of the Friedel oscillation, determines the value of the 
SDM. That is why there is an excellent agreement between SC-SJM and FPLAPW work 
functions for large n, in spite of the fact that our Al(lOO) slabs contract while Al(lOO) 
slabs expand in FPLAPW. To explain the differences between our results and those of 
FPLAPW in the work function of Al(llO) slabs, we notice from the bottom sub-figure 
of Fig. Hlthat FPLAPW predicts more compression than SC-SJM for n=2, and hence, 
smaller Fermi wavelength which leads to a rapid bending of the work function curve. For 
n > 4, we notice a constant phase difference between the two oscillatory curves. That 
is, (aside from the vertical shift which we discussed) the value of the SDM at n=4 in 
SC-SJM is equivalent to that in FPLAPW at n=6, and so on. This fact can be explained 
by assuming that the overlap between the major peaks in the Friedel oscillation at both 
sides of the slab effectively vanishes for n > 6 atomic layers in the FPLAPW and for 
n > 4 in the SC-SJM. For larger n's, the constant phase difference means that the Ap 
of the electrons in the bulk region is the same in the SC-SJM and the FPLAPW. Now, 
we focus on the top sub-figure of Fig. [6] on Al(lll). For n > 7, aside from the vertical 
shift, the two curves have the same frequency of oscillations and therefore, in the bulk 
region both SC-SJM and FPLAPW have the same electronic Ap. However, for n < 7, 
from n=2 to n=3 FPLAPW gives contraction while from n=3 to n=7 it gives expansion 
(See top sub- figure in Fig. H]). Therefore, going from n=2 to n=7, first Ap decreases 
and then increases. As a result, in the FPLAPW work function plot we observe a rapid 
bending followed by a slow bending. From the FPLAPW results, it can be deduced 
that at n=7, the surface formation is completed. To sum up, (from Fig. Hj) FPLAPW 
results show positive slopes at n=2, followed by a negative one for Al(lOO) and Al(llO), 
whereas it shows negative slope at n=2 followed by a positive one for Al(lll). These 
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behaviors are reflected in the work function curves (Fig. E]) as slow bending at n=2 
for Al(lOO), Al(llO), and rapid bending for Al(lll). The dotted straight lines in each 
sub-figures correspond to the result obtained from the self-consistent SJM calculations 
for semi-infinite system [T9l [20] . For completeness, we have also included the PP-LDA 
results[2l[Hl[7]. 

To see the effects of surface corrugation on the work functions, we have used the 
face-dependent relation of Eq. ([8]) for {6v) for a selected Al(llO) slabs. The results are 
presented in Table El 



Table 2. Comparison of the relaxations, averaged difference potentials, and work 
functions of n-layered Al slabs as well as semi-infinite jellium results with those of 
FPLAPW. 



n 


case 


AL/L{%)^ 


(Sv) (eV)^ 


W (cV)^ 


AL/L(%)b 


{Sv) (eV)^ 


W (eV)^ 


AL/L{%y 


W (eV)'= 


9 


(100) 


-1.12 


-2.45 


4.18 


-2.67 


0.20 


3.69 


2.1 


4.23 




(110) 


-2.19 


-2.42 


4.13 


-6.55 


3.13 


3.56 


-2.4 


4.16 




(111) 


-0.91 


-2.46 


4.27 


-1.24 


-1.70 


4.09 


3.6 


4.05 


11 


(100) 
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-2.46 


4.26 


-2.02 


0.18 


3.78 


5.7 


4.28 




(110) 


-1.41 


-2.44 


4.27 


-5.32 


3.09 


3.61 


-12.6 


4.00 




(111) 


-0.58 


-2.47 


4.25 


-0.86 


-1.71 


4.07 


4.4 


4.07 


13 


(100) 


-0.56 


-2.47 


4.27 


-1.56 


0.17 


3.77 


7.1 


4.28 




(110) 


-1.13 


-2.45 


4.26 


-4.44 


3.06 


3.62 


-2.3 


4.09 




(111) 


-0.38 


-2.47 


4.24 


-0.61 


-1.72 


4.06 


1.5 


4.03 


15 


(100) 


-0.40 


-2.47 


4.24 


-1.25 


0.16 


3.74 


2.7 


4.25 




(110) 


-0.81 


-2.46 


4.23 


-3.64 


3.03 


3.54 


-7.4 


4.06 




(111) 


-2.08 


-2.48 


4.26 


-0.41 


-1.73 


4.08 


3.3 


4.06 


20 


(100) 


-0.09 


-2.48 


4.26 


-0.07 


0.14 


3.76 








(110) 


-0.44 


-2.47 


4.24 


-2.44 


3.00 


3.55 








(111) 


-0.02 


-2.48 


4.24 


-0.01 


-1.74 


4.05 






semi-inf^ 


flat 
(100) 




-2.49 


4.24 




0.12 


3.62 








(110) 










2.92 


3.81 








(111) 










-1.74 


3.72 







^ Present work, flat surface. 
^ Present work, face-dependent corrugation. 
FPLAPW (Ref [TD]). 

Semi-infinite SJM system from Ref. [19]. This result is quite close to that in Ref. ^20] which is 4.27 
eV. 



In Table the numbers in the fifth column are the work functions of flat surface 
slabs. The eighth and tenth columns are the face- dependent and the FPLAPW results, 
respectively. As is observed from the Table [21 the surface corrugation term has affected 
the amount of relaxation and the averaged potential difference, (Sv). Any change in 
the difference potential, causes a change in the SDM, and thereby lead to a change 
in the work function. However, the (100) work functions of fiat surface agrees nicely 
with those of FPLAPW[10j, and therefore, the face- dependent term worsens the results. 
On the other hand, work functions of both (110) and (111) slabs with fiat surfaces are 
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20 



Figure 7. Same as in Fig. [6] with semi-infinite results 2.92 eVpTO which is close to 
2.94 eVEni. 



higher than those of FPLAPW, but, the corrugation term causes overcorrection and 
worsening the agreement with FPLAPW. The lowest four rows of the table show the 
results obtained from application of the SJM to semi-infinite jellium system [T5] which 
are in excellent agreement with our results for large n's. 

In Fig. [7] we have presented the work functions of Na slabs at their equilibrium 
states for (110), (100), and (111) surfaces as functions of number of layers. The dotted 
lines correspond to the SJM results (2.92 eV) for semi-infinite system[Tn] which is close 
to 2.94 eV obtained in Ref. [20]. 

The calculation results for the work functions of Cs slabs are plotted in Fig. [81 The 
dotted lines correspond to the SJM results (2.24 eV) for semi-infinite system [T9] which 
is close to 2.26 eV obtained in Ref. [20j. As in the case of rj results, here also, the 
behaviors are similar to those of corresponding Na slabs. 

In Fig. [9l we have plotted the total energies per electron for Al(lll), Al(lOO), and 
Al(llO) slabs at their equilibrium states. As seen, the behaviors are quite smooth and 
approach to the bulk value -19.1 eV [ilSj from above for large enough L's. Similar 
behaviors are obtained for Na and Cs slabs with respective asymptotic values of -6.26 
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Figure 8. Same as in Fig. [6] with semi-infinite results 2.24 eVfT9' which is close to 
2.26 eV^. 



and -4.64 eV 

The surface energy of the slab is defined by 

^iL^) = ^[Esm{L^)-Elf^{L)], (10) 

where EsmiL^) and ^sjm(^) are the SC-SJM and bulk-SJM total energies of the slab 
with surface area A. The bulk-SJM total energy is proportional to L = ndhki and 
Eq. (fTOj) can be rewritten in the form of 

where n and e^"''^ are the bulk density and bulk SJM energy per particle, respectively. 
We have fitted the values of the SC-SJM total energy per unit area to Eq. (ITT!) to obtain 
^buik^ and used it in Eq. ffTOj) to obtain cr(L^). In the fitting we have used the results 
with n > 6 to eliminate the significant fiuctuations and obtain more accurate values for 

^bulk 

In Fig. [TOl we have plotted the surface energies obtained from using Eq. (ITOl) for 
Al slabs with fiat surfaces. The results are compared with those obtained from the 
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Figure 9. Total energies per electron of Al slabs at their equilibrium states for (111), 
(100), and (110) surfaces. The dotted lines correspond to the bulk value, -19.1 eV [T3]. 



FPLAPW[TO]. Comparison shows a systematic underestimation of ~ 8% for Al(lll), 
~ 9% for Al(lOO), and ~ 10% for Al(llO) surfaces, respectively, which can be partly due 
to the GGA used in the FPLAPW calculations. In the bottom sub-figure, Al(llO), we 
have also included the PP-LDA results which lie above both the SC-SJM and FPLAPW 
results. As is observed, the structures of the plots are in excellent agreement with those 
of FPLAPW. On the other hand, it is seen that both PP and FPLAPW plots have the 
same structures. The disagreements of n=2 in Al(lOO) and in Al(llO) with those in 
FPLAPW are consistent with those in Fig. [61 The phase difference is also observed in 
the bottom sub-figure. The dotted and dashed lines correspond to the results of SJM 
for semi-infinite system with 953 erg/cm^ |19] and 925 erg/cm^ |20j. respectively. 

The surface energies of Na and Cs slabs are plotted in Figs. [11] and [121 respectively. 
Consistent with the top sub-figures in Figs. [21 [3], [TJ [HI the top sub-figures here also show 
smooth behaviors. Compared to any FPLAPW calculations, we expect a systematic 
underestimation of surface energies of at most ~ 10% for Na and Cs as well. On the 
other hand, as in Al case, our results for Na and Cs are in excellent agreement with 
those of the semi-infinite SJM for fiat surfaces [T9l [20] . Finally, we note that the semi- 
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Figure 10. Surface energies of Al slabs at tlieir equilibrium states for (111), (100), 
and (110) surfaces. The labels are the same as in Fig.[Sl The dotted and dashed lines 
correspond to the SJM results for the semi- infinite system with 953 erg/cm^ T9] and 
925 erg/cm^ [501, respectively. 

infinite calculations [13] with surface corrugations have raised the Al surface energies to 
977, 1103, and 921 erg/cm^ for (100), (110), and (111) indices, respectively, which are 
in better agreements with those of FPLAPW. 

4. Conclusions 

We have applied the SC-SJM to predict the electronic properties of simple metal slabs of 
Al, Na, and Cs at their equilibrium states. To make a direct correspondence to atomic 
slabs, we have considered only those L values that correspond to n-layered atomic slabs 
with 2 < n < 20 for surface indices (100), (110), and (111). In our calculations, the 
slabs are taken as isolated systems which can undergo relaxations in the z direction. 
Our results show that, in contrast to the r* values (obtained in Ref. |9]) which are 
always smaller than those of the bulk, for some slabs the r]'s assume values that are 
larger than those of the bulk. That is, in some cases the relaxations are realized as 
expansions. Our results show that all Al slabs are contracted while some of Na and Cs 
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Figure 11. Same as in Fig. [10] for Na with 171 erg/cm^ for semi-infinite SJM[19l[20] 



slabs expand. From these results we explore that if the interlayer spacing happens to be 
so small or the Fermi wavelength happens to be so large that the condition 2Af /d^ 1 is 
met, then there would be a chance that the QSE be realized as consecutive expansions 
and contractions. Since this condition is not met in the Al case, we do not see any 
significant sharp oscillations in the QSEs. For more justification of the criterion, we 
have tested on Li slabs for 2 < n < Q and have obtained some expansions. This 
criterion along with the first-principles or experimental results can be used to estimate 
an effective free-electron value for a given slab. On the other hand, our results for the 
work functions and surface energies show similar structures as in those of the FPLAPW 
results. We have explained the quantitative differences between our results and those of 
the FPLAPW to be partly due to the fact that our model calculations overestimate the 
SDM values, and partly because of the fact that in FPLAPW calculations, the GGA 
has been used which give significant improvements over the LDA results for systems 
with high inhomogeneities in electron densities. Finally, our results are in excellent 
agreement with the self-consistent results of the SJM for semi-infinite system. 
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Appendix A. Relation between rj and r* 

The energy of the slab system is a function of volume, V, and r^, i. e., E = E{V,rs). 
In fact, instead of V, it depends on L because, in the slab system the surface area 
does not change. However, we use V for generality of our arguments. We define two 
different derivatives, Dn and Dy, which are derivatives when is constant and when 
V is constant, respectively. That is. 



D - ^ 
dr.^ 



E{V,rs 
N 



N 



1 d 

Ndr\ 



Att d 

-rJ^|B(V.r,)]„. (A.1) 
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d 

dr^ 



E{V,r, 
N 



1 d 
Vd^, 



V_ 

N 



E{V,rs) 



V 



An d 
~ 3Vd^s 

However, on the one hand, using Eq. ([7]) one obtains 

dV _ 3V 

drs rs ' 
which results in 



fdE\ dV 

[dvj^^d^,^ 



dE\ 
drs J 



N 



dV 



drsj , 



and on the other hand, the constant volume conditions results in 



Now, combining Eqs. (lA.ip - (]A.5p one obtains 



Dv{V,rs) = D^{V,rs) + 



Anrl 



E{V,rs 



V 



fdE_ 

\dV 



By definition. 



Dy(\/,r:) = 0, 



and 



DN{y,rl)=Q. 
Inserting Eqs. (IA.7p and flA.Sp into flA.6p results in 



DN{V,rl) + 



V 



and 



V 



0. 



(A.2) 



(A.3) 



(A.4) 



(A.5) 

(A.6) 

(A.7) 
(A.8) 

(A.9) 
(A.IO) 



Taylor expansion of Eq. flA.QP around r\ to linear term and using Eq. flA.SP and 
combining with flA.lOp one obtains 

^ . Dv{r\) 



= 



dPNirl) dDvjrl) 
drl drl 



(A.ll) 



Since the energy E/N is a convex function of r^, both terms in the denominator of 
Eq. (lA.lip are positive while the sign of the numerator is positive for r* < r] and 
negative for r* > r]. 
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